Electron-hole imbalance in superconductor-normal metal mesoscopic structures. 
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We analysed the electron-hole or, in another words, branch imbalance (BI) and the related elec- 
tric potential Vimt which may arise in a mesoscopic superconductor/normal metal (S/N) structure 
under non-equilibrium conditions in the presence of a supercurrent. Non-equilibrium conditions 
can be created in different ways: a) a quasiparticle current flowing between the N reservoirs; b) a 

■ temperature gradient between the N reservoirs and no quasiparticle current. It is shown that the 
voltage Vimb oscillates with the phase difference Lp. In a cross-geometry structure the voltage Vimb 

(— I arises in the vertical branch and affects the conditions for a transition into the tt— state. 
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^ ' A few decades ago a great deal of interest was paid to the study of effects related to the so called branch imbalance 
I (BI) (see Refs. [Q, || and for example the reviews [|[ The BI implies that populations of the electron- like and 
hole-like branches of the excitation spectrum in a superconductor or a normal metal are different. For example, the 
BI may arise in a superconductor near the S /N interface if a current flows through this interface and the temperature 
is close to Tf.. The conversion of the quasiparticle current Jq into the condensate current js occurs over a rather 
, long length Aq called a BI relaxation length. Over this length populations of the electron-like and hole-like branches 
of the energy spectrum differ from each other. The difference between these populations is characterized by the 
O I distribution function /_ = — (n| —pi); the function is the distribution function of the electron-like excitations and 
Pl{e) = 1 — n|(— e) is the distribution function of the hole-like excitations. In the considered case of a spin-independent 
interaction one has n| = and p| — Pi- One can show that the function /_ differs from zero if divjq 7^ 0. The 
non-zero distribution function /_ leads to the appearance of an electric potential Vimb in a superconductor (or in 
a normal metal) which can be expressed in terms of the function /_ (see below). The BI may also arise in a bulk 
qI^ ■ superconductor. For example, if longitudinal collective oscillations with a finite wave vector q are excited in the 
^) ' superconductor, the BI arises because in this case divjq = iqjq ^ . When these modes are excited (they are weakly 

■ damped only near T^), the quasiparticle current jq oscillates in a counter phase with the condensate current jg, so 
' that the total current remains equal to zero. These oscillations have been observed experimentally by Carlson and 

Goldman (Carlson-Goldman mode) ^ and have been explained theoretically in Refs. ||, 0. Another example of a 
■»«^ ] system, in which the BI arises, is a uniform superconducting film in the presence of a temperature gradient VT and 
"j^ ■ a condensate flow. It was established experimentally Q and theoretically |lO|, |ll| that in this case the BI has a 
' magnitude which is proportional to WsVT, where Vs is the condensate velocity. 
I I Recently there has been growing interest in the study of transport properties of S/N mesoscopic structures. Several 
' O \ interesting, phase-coherent effects have been observed in these structures. Among them one can mention the change 
^ of sign of the Josephson critical current Ic in a four-terminal S/N/S mesoscopic structure. If an additional dissipative 
current (or an applied voltage) between the N reservoirs in a S/N/S structure of a cross geometry exceeds a certain 
value, the current Ic changes sign and the Josephson junction is converted into the 7r-state( a theory for this 
effect was developed in Refs. ||,]l5|). 

In this Letter, we consider S/N mesoscopic structures under non-equilibrium conditions and study the BI which can 
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H \ arise in the N conductor if there is a flnite phase difference tp between the superconductors in these structures. We 
. . 1 briefly analyse effects arising due to the BI and discuss how some of them can be observed. To our knowledge these 
effects are missed in most papers on transport properties of mesoscopic S/N structures. In some cases the physics 
of the BI in mesoscopic S/N structures resembles that of the BI in a uniform superconductor under nonequilibrium 
conditions. However even in these cases the BI has its own specific characteristics. For example, the voltage Vimbi 
which is associated with the BI, is an oscillating function of the phase difference between the superconductors and 
it may appear even in the absence of a temperature gradient. Consider for example the structure in fig. |l|a. We will 
show that the voltage Vimbiy) (or the electric field) arises in the vertical wire regardless of how the system is brought 
out from the equilibrium state. One can apply a temperature gradient between the N reservoirs disconnected from 
the external circuit or one can pass a finite current between the N reservoirs as it was done in the experiment [E2[. 
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FIG. 1: Schematic view of the structures under consideration. Barriers or bad S/N interfaces are represented by thick lines. 
The dashed and dotted lines show the spatial distribution of the applied voltage (between the N reservoirs,fig. |l^) and the 
voltage caused by the charge imbalance. The electric potential in the superconductors is assumed to be zero. 



When considering S/N/S structures, it is convenient to introduce two types of the distribution functions: f± as it was 
done in Ref. The function /+ determines the order parameter in the superconductors and, for example, the 

condensate current in S/N/S structures. This function is related to the distribution functions of electrons n| = 
and holes pi — pi: /+ = 1 — (n| + Pi). The Josephson current in the N wire can be written in the form 



Js 



ia/2)W / def+.L 



(1) 



where W is the cross-section area of the y-wire, the distribution function /+, generally speaking, deviates from its 
equilibrium form /eg = tanhe/3 and should be determined from the kinetic equation, here /3 = 1/2T. The function 
Js is "a partial condensate current" and is expressed in terms of the condensate Green's functions F^(^): = 
{l/4)Tr{a^{F^dyF^ - F^dyF^)}. This "current" does not depend on y. On the other hand, as follows from the 
conservation of the current, js can be written as a current through the S/N interface. This means in particular that 
"the current" Js can be represented in the form Js = {rg/ Ly)Tr{az{F^F^ — F'^Ff )}, where the functions F^*^'^) 
must be taken at the S/N interfaces; = Ry/Rs, Ry — gLy is the resistance of the y-branch (per unit area) and Rs 
is the the resistance of the S/N interface (per unit area). Under nonequlibrium conditions, which may have a different 
origin (non-equal temperatures of the S and N reservoirs, an additional current between the N reservoirs etc), the 
current given by eq. (|l|) is not equal to the Josephson current through the S /N interface because the last one depends 
not only on the value of at the interface but also on the distribution function fs+ in the S reservoirs which is 
assumed to be equilibrium (/5_|_ — feq). This means that an electric field E = —dyVimb arises in the y-branch even 
if there are no voltage differences neither between N reservoirs nor between the S reservoirs. The electric field drives 
the quasiparticle current which compensates the above mentioned difference between the condensate currents. The 
appearance of the electric field and the quasiparticle current means that strictly speaking eq. (|l]) does not describe the 
maximum current in the y-branch in the absence of the voltage difference between the S reservoirs (in most papers on 
this subject the maximum current was found from eq. (|l|)). In order to find the critical current in a non-equilibrium 
situation, one has to use a more general formula. However at low temperatures eq. (^ determines the critical current 
Ic with a good (exponential) accuracy. 

In the present paper we consider S/N mesoscopic structures of two configurations (see fig. |^). In simple limiting 
cases we find the distribution functions and analyse the BI arising under non-equilibrium conditions. The distribution 
functions f± obey the kinetic equations |16 which for the structure shown in fig. |l|b can be written in the form (see 
for example p^) 



Ld4M± d^f±{x) +Jshix) ± Jan d^f^ix)] =rs[A± Six - L, )+ A± Six + L, 



(2) 



where all the coefficients are expressed in terms of the retarded (advanced) Green's functions: G^''^-' = G-^^'^^^^'ctz + 
pR{A). ^ (1 _ qRqA ^ (i?fli?A)^)/2; Jan = {F^F^)z/2, Js = (l/2)(F«9,i?« - F^dJ^)z, A± = (yvs + 

5iT)(/±-/s±)-(ff.±/sT+ff.T/T); 5i± = (l/4)[(F«±F^)(F|±F#)]i; = {Xj 4)\{F^ TF^){Ff ±Ff%. The 
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parameter rg = Ri/Rs is the ratio of the resistance of the N wire Ri = pLi and S/N interface resistance Rs; the 
functions v4_ and A+ coincide with A-,A+ if we make a substitution tp —Lp. We introduced above the following 
notations {F^F'^)i = Tr{F'^F'^)/2, (F^F^)^ = Tr{a^F^F'^)/2 etc.; v, vs are the dcnsity-of states in the N film at 
X — Li and in the superconductors. The functions fs± are the distribution functions in the superconductors which 
are assumed to have the equilibrium forms. This means that fs+ = feq ~ tanh(e/3o) and fs- — 0, because we set the 
potential of the superconductors equal to zero. We neglect branch imbalance in the superconductors assuming that 
the distribution functions f± recover quickly their equilibrium forms in S due to a big size of the superconductors 
in comparision with the size of the S/N interface. In the case of the structure in fig. |l|a the left-hand side of eq. (|^) 
should be written down for each branch of the structure and be set equal to zero. At the S/N and N/N' interfaces 
we use the boundary conditions which are given by the right-hand side of eq. (H) (at the N/N' interfaces the index 
S should be replaced by the index N' and all the condensate Green's functions in N' should be set equal to zero). 
Consider the structure shown in fig. |^a. Eqs. (|l|) can be integrated once in each branch, and, for example, in the 
y-branch we obtain 

M± dyf±{y) +Jsh{v) ± Jan dyf^iy) = Jy^. (3) 

where Jy^ are the total "partial currents" in the upper (y > 0) and lower {y < 0) parts of the j/-branch of the N 
wire. The current Jy^ is a "vector" with the components {Jy+, Jy- )■ At the crossing point the current conservation 
law takes place Jj,> -I- Jx> — Jy< + Jx<- At the S/N interface the current Jj^> is related to the Green's functions at 
y = Lv 

Jy> = {rs/Ly) A. (4) 

In order to obtain the current through the y-wire, the "current" Jy> should be substituted into the integrand of 
eq. (|^ instead of f+Js- We present the solutions for the distribution functions assuming first the weak proximity 
effect. This means that the amplitude of the condensate functions in the N wire should be small: << 1. In 

this case F^^^^ can be easily found from the linearized Usadel equation (the solution for the structure shown in fig. ^ 
is presented in ref. As follows from the form of the functions F^''^\ they are small if the condition e >> Syrs 

is satisfied, where the characteristic energy s is equal to the Thouless energy Sy = D/Ly ( in the case of the geometry 
in fig. ^ the Thouless energy is £i = D /L\) or to the temperature T. In the case of the weak proximity effect, a 
solution for the kinetic equations (□) also can be easily found with the help of expansion in the parameter rs (for 



the case > 1 eq. (El) was solved numerically in Ref. |14 ). We consider two types of the nonequalibrium situation: 



a) the temperatures of all the reservoirs are the same, but the electric potentials at the left and right N reservoirs are 
±y (a current flows between these reservoirs); b) no current between the N reservoirs, but the temperatures of the 
left and right N reservoirs are different, so that the distribution function in the right (left) reservoirs F^^^is equal to: 
Fr^i = tanh(e/3r.i)- I^i the main approximation we find 

a)f+{x) - Fv+-J+{v)^Fv+ 

h)f+{x) = {Fi+Fr)l2 + {x/Lx){Fr-Fi)/2-J+{y) = {Fi+Fr)/2 (5) 

and 

a)f-{x) = {x/Lx)Fv-;f-{y) ^rs{y/Ly)g2^{feq - Fv+); 

&)/_(x) = 0;/_(y) = r5(2//L,)5,_[/e, - {Fr + Fi)/2] (6) 

Knowing the distribution function /_, we can calculate the electric potential in the y-branch V with the help of 
the formula M, Il6| 



eViy) - (l/8)Tr J deG{e, y) = (1/2) j dev{e)f^ (7) 

where G = G^f — fG^ is the Keldysh component of the matrix Green's function, / = 1/+ -t- a^f-. is the matrix 
distribution function, u(e) = {l/A)Tr(az{G^ — G^)) is the density-of-states in the N wire. We easily find from eqs. (^ 
1^ that the potential V{y) is an odd function of y and the electric field E{y) is an even function of y. We also see that 
the electric field arises only if the phase difference between the superconductors is not zero; otherwise the function 
gz- is zero. The field E or the potential V{y) oscillate with increasing phase difference if. As is seen from eqs. (^0), 
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FIG. 2: The dependence of the critical current 7c (in arbitrary units) on the applied voltage V for two temperatures: T — 0.5eL 
and T — 2e_L, where e_L — D/L^ ,Ly = 0.3L. The value e_t is chosen in such a way that A(0) = 5eL- The solid lines are calculated 
from the correct eq. (^) and the dashed lines are calculated from the approximate, widely used eq. (^) . The difference between 
these curves determines the voltage Vy (divided by the S/N interface resistance) at the S/N interface. 



the electric field arises regardless of the origin of the non-equiUbrium state: the function /+ may deviate from the 
equilibrium function feq if a finite current flows between the N reservoirs or if the temperature of the N reservoirs 
differs from the temperature of the S reservoirs, fn the considered case of a weak proximity effect the conductance G 
between the N reservoirs decreases with increasing the phase difference ip (at small ip). We also have calculated the 
critical (maximum) current using the correct expression eq. (^) and the approximate one eq. (0) (see fig. |^) . 

We see that the difference between two curves is significant if the temperature is not low. The difference between 
the critical currents is determined by the quasiparticle current jg = aE. This follows directly from eq. (^ if one takes 
into account that in the main approximation A/_ = 1 and the last term on the right-hand side can be neglected, ft is 
worth noting that /_ 7^ only in the normal wire. In the superconductors the function /_ and therefore the voltage 
V are assumed to be zero, so that at the S/N interface there is a voltage drop from a finite value of V in the N wire 
to zero in the superconductor. 

We also considered another limiting case when the proximity effect is not weak, that is, the condensate function in 
the N wire |F| is not small. The obtained results qualitatively are similar to those which have been established for 
the weak proximity effect. The only difference is that the conductance in this case increases with increasing phase 
difference (p. Different behaviour of the conductance G as a function of ip was studied in Refs. where a transition 
from a decreasing to increasing dependence (at small (p) of G{ip) was obtained by varying the applied voltage. In a 
recent paper [|l9j a similar transition (obtained by varying the temperature) was studied in detail both experimentally 
and theoretically. Contrary to Refs. it was assumed in Ref. that the S/N interface is perfectly transparent. 
A good agreement between theoretical results and experimental data was obtained. 

At last we consider a specific thermoelectric effect arising in the structure in fig. |l|b. This effect was measured 
recently in mesoscopic S/N structures As established in j23|, if temperatures of the normal reservoirs 

are different, a voltage arises between these and superconducting reservoirs (the S reservoirs have the same electric 
potential because they are connected with a superconducting loop). The origin of this voltage which can be called 
thermoemf is completely different from the ordinary thermoemf in S/N/S junctions studied in Refs. p^ . In the 
last case the thermoemf appears due to the ordinary thermoelectric component of the quasiparticle current which is 
neglected here. In the limit of the weak proximity effect we have calculated the distribution functions f± and the 
voltages Vi^r caused by the temperature difference between the N reservoirs; here Vi^r are the electric potentials at 
the left and right N reservoirs, respectively (the electric potential at the S reservoirs is set to zero). We assumed that 
the temperature of the left N reservoir To coincides with the temperature of the S reservoirs and the temperature of 
the right N reservoir T is elevated: T — To+ 6T. The distribution functions can be easily found using an expansion 
in the parameter rs- For the voltages V± = {Vr ± V/)/2, we obtain from eq. (0) 



eV+ 



eV. = rs5T{L^/2L) J de(e/3)ff,_(e, Li)^ (8) 
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FIG. 3: The temperature dependence of the vohages V+ = {Vr + Vi)/2 and V- = {Vr — Vi) caused by an extraordinary, 
phase-coherent thermoemf for different values of the parameter A(0)/ez,: A(0)/ez, — 1.0; f .4; f .8; 2.2; 2.6; 3.0; 3.4 and 4.0 for the 
curves 1 — 8, respectively. 



where /^^ = cosh~^(e/?). The expression for V+ was obtained in Ref. ||2^. Here we also obtained the voltage difference 
between the N reservoirs V-. One can see that this voltage contains the small parameter rs in comparision to 
the voltage V+. Both voltages are proportinal to sin ip, that is, they oscillate with increasing phase difference. The 
temperature dependence of the amplitudes of V± are plotted in fig. ^. 

We see that the maximum of the amplitude of V+ is located at a lower temperature than the maximum of V- . In 
principle this behaviour may lead to a nonmonotonic behaviour of the voltages Vr or VJ as a function of temperature. 
The change of the phase of the voltage oscillations with increasing phase ip may have the same origin as that in the 
case of the conductance i.e. the change of the dependence ip{H) , where H is the external magnetic field. Our 

results can not be compared quantitatively with the recent experimental data [2l| | because the experimental structure 
corresponds to the case ^ 1. We will analyse this more complicated case in a separate paper. 

In conclusion, we have analyzed the branch imbalance effects in S/N mesoscopic structures. We have shown that in 
the structure in fig. [^a a voltage Vy related to the BI is set up if a current is driven through the x-axis or a temperature 
gradient exists between the N reservoirs. The voltage Vy is proportional to sin 95, i.e. it oscillates with increasing 
phase difference ip. A similar voltage Vimb arises in the vertical wire if a temperature gradient ST exists between the N 
reservoirs. We also studied an unusual thermoelectric effect in the structure shown in fig. In this case voltages 
arise in the right and left N reservoirs (the electric potential at the superconductors is assumed to be zero) if there is 
a temperature difference between the normal resevoirs. These voltages are proportional to ST simp and they are not 
related to the ordinary thermoemf because we have ignored the small thermoelectric component in the quasiparticle 
current. 
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